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ABSTRACT 


An  experimental  and  analytical  program  was  conducted  as  a  step  in  a  syste¬ 
matic  engineering  parametric  study  to  establish  the  necessary  engineering 
criteria  for  the  prediction  and  prevention  of  flow  decay  in  operational 
nitrogen  tetroxide  systems.  Tests  were  made  with  both  green  and  red-brown 
nitrogen  tetroxide  over  ranges  of  seven  other  parameters:  initial  propel¬ 
lant  temperature,  temperature  drop  before  reaching  the  test  section,  time 
during  which  the  temperature  drop  is  imposed,  filter  pore  size,  local  velo¬ 
city  through  the  filter,  total  volume  per  unit  filter  area  of  propellant 
passing  through  the  filter  during  a  test,  and  iron  saturation  condition. 

Flow  decay  was  observed  under  some  test  conditions  and  was  absent  under 
other  conditions.  The  rates  of  decay  ranged  up  to  2.8  %/min,  although  the 
mean  rate  was  about  0.5  ?o/min  for  all  tests  in  which  any  significant  flow 
decay  occurred. 

It  was  possible  to  deduce  information  about  the  effects  of  these  eight  major 
independent  variables  on  flow  decay.  However,  it  was  found  that  the  effects 
of  these  parameters  are  generally  not  simple  or  independent  of  each  ether; 
they  exhibit  many  interactions.  In  addition,  there  are  often  threshold  ef¬ 
fects  for  flow  decay  (i.e.,  an  identifiable  boundary  between  a  range  of 
variables  for  which  no  flow  decay  occurs  and  a  range  for  which  flow  decay 
occurs  at  varying  rates).  These  thresholds  are  not  sharp,  and  further  de¬ 
pend  upon  the  interactions  of  the  independent  variables. 
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NOMENCLATURE 


A  Filter  area 

A^  Free  area  of  filter  available  for  flow 

C  Discharge  coefficient  (defined  by  Eq.  8) 

D  Filter  nominal  pore  size 

H  Pressure  drop  across  the  filter,  psi 

H*  Calibration  value  of  H  (i.e.,  for  a  clean  filter),  psi 

k  Constant  in  Eq.  3 

L  Flow  decay  parameter  defined  by  Eq.  11,  gal./sq  in. 

Flow  decay  parameter  defined  by  Eq.  10,  gal./sq  in. 
n  Constant  in  Eq.  3 

Q  Flowrate,  gal. /min 

Q*  Calibration  flowrate  (i.e.,  for  a  clean  filter),  gal. /min 

Q2  Best  estimate  of  the  value  of  Q2  (i.e.,  flowrate  at  time  tp  if 

H2  *  H1 

Re  Reynolds  Number  fcr  filter,  DQp/Ay 

t  Time,  min 

T  Temperature,  F 

Tq  Initial  propellant  temperature,  F 

AT  Temperature  drop  (T  minus  temperature  at  test  section  during 

flow),  F 

V  V^/A  (total  amount  of  propellant  per  unit  filter  area  which  has 
passed  through  the  filter),  gal./sq  in. 

Volume  of  flowing  through  filter  between  t^  and  t7,  gal. 

W  1  -  Q/Q*,  fractional  drop  in  flowrate  between  the  clean  filter 

case  and  the  case  of  interest 
X  1  -  Q^/Qp  fractional  flow  decay  at  constant  H 

y  General  random  variable 

Y  X/(t2  -  t,) 

z  (w2  -  wx)/(t2  -  tx) 


GREEK 


e 

u 

u 

p 

0 


Random  error  term 
Viscosity 

Mean  value  of  the  i^  variable;  true  value  of  the  outcome  of  the 
.th 

1  experiment 
Density 

Standard  deviation 


LITERAL 


cov(a.b) 

E(y) 

GB 

RL 

SAT 


var (a) 


Covariance  of  a  and  b 
Expected  value  of  y 

Two-level  variable  with  assigned  values  of  1  when  green  is 

used  and  -1  when  red-brown  N-0.  is  used 

2  4 

Two-level  variable  with  assigned  values  of  1  for  rapid  cooldown 
and  -1  for  long  cooldown 

Two-level  variable  with  assigned  values  of  J.  for  N204  which  is 
doped  with  iron  pentacarbonyl  to  ensure  iron  saturation,  and  -1 
for  N',,0.  tested  as  received 
Variance  of  a 


SUBSCRIPTS 


o 

1 

2 


Initial  value 
Value  at  time 
Value  at  time 


at  t  =  o  as  with  clean  filter 
C1 

t2  (t2  >  .,) 
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INTRODUCTION 


The  phenomenon  of  flow  decay  is  defined  as  a  spontaneous  decrease  in  the  flow- 
rate  through  a  constant  pressure  flow  system.  Flow  decay  in  nitrogen  tetrox- 
ide  systems  is  a  function  of  the  in  situ  formation  of  solid  or  gel-like  mate¬ 
rials  that  can  obstruct  the  flow  through  valves,  filters,  orifices,  or  any 
other  flow  element  with  a  constriction  of  small  size.  There  are  several 
ferric  nitrate  derivatives  which  can  be  deposited  to  produce  flow  decay,  such 
as  nitrosyl  tetrani  tratoferrate ,  NOFeCNO-j)^,  and  partially  hydrated  or  hy¬ 
drolyzed  ferric  nitrates,  Fe(N03)3_n(0H)n  .  xlHO.  These  materials  are  solu¬ 
ble  in  nitrogen  tetroxide  at  levels  of  the  order  of  a  few  parts  per  million 
(as  iron).  Nhile  there  are  many  factors  which  can  influence  the  corrosion 
processes  that  produce  these  compounds  in  solution,  the  extensive  exposure 
of  all  nitrogen  tetroxide  to  iron  during  the  manufacturing  and  shipping  proc¬ 
esses  can  be  expected  to  form  these  materials  in  approximately  equilibrium 
amounts . 

Rocketdyne  has  conducted  a  number  of  previous  investigations  of  aspects  of 
flow  decay,  beginning  in  1964  (Ref.  1  through  4).  Laboratory  investigations 
and  observations  of  flow  decay  in  field  systems  have  also  been  made  by  a  num¬ 
ber  of  other  companies  and  agencies.  A  useful  summary  of  past  work  is  given 
in  Ref.  5.  The  previous  studies  have  provided  some  understanding  of  the 
underlying  mechanisms  and  many  of  the  important  parameters  affecting  nitrogen 
tetroxide  flow  decay.  However,  in  general,  they  have  not  provided  informa¬ 
tion  that  could  permit  prediction  of  the  occurrence  or  absence  of  flow  decay. 
The  fulfillment  of  this  technology  void  would  provide  the  basis  for  engineer¬ 
ing  control  of  this  potential  system  failure  mode  through  Jciimtion  of  re¬ 
quired  system  design  criteria  and/or  system  management  concepts. 

This  program  was  intended  to  be  the  initial  step  in  a  systematic  engineering 
parametric  study  to  establish  the  necessary  engineering  criteria  for  the 
prediction  and  prevention  of  flow  decay  in  operational  nitrogen  tetroxide 
systems.  Its  goals  were  to  establish  experimental  techniques  ar.d  data  analy¬ 
sis  techniques  to  give  statistically  reproducible  nitrogen  tetroxide  flow 
decay  data,  and  to  provide  an  interim  guide  to  the  gross  effects  of  major 
parameters  on  flow  decay.  A  fundamental  assumption  underlying  the  pro¬ 

gram  was  that  any  commercial  nitrogen  tetroxide  will  have  the  potential  of 
undergoing  flow  decay,  if  suitable  conditions  are  present  during  flow. 
Therefore,  the  philosophy  was  to  disengage  corrosion  studies  (i.e.,  varia¬ 
bles  which  primarily  affect  the  formation  of  iron  compounds  in  solution)  from 
these  investigations  of  deposition  of  flew  decay  material. 

In  this  program,  the  first  task  (Phase  I)  was  to  develop  a  test  plan  to  be 
used  in  performing  flow  decay  tests.  In  the  second  task  (Phase  II),  flow 
techniques  and  apparatus  developed  and  used  in  previous  Rocketdyne  programs 
were  modified  and  used  to  perform  parametric  tests,  following  the  sequential 
experimental  design  outlined  in  Phase  I.  The  third  task  (Phase  III)  comprised 
the  analysis  of  test  results  to  define  sets  of  conditions  under  which  flow- 
decay  will  or  will  not  occur  in  specific  system  components,  and  to  develop  a 
preliminary  basis  for  establishing  criteria  to  avoid  this  failure  mode. 
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FLOW  DECAY  VESTS 


EXPERIMENTAL  SYSTEM 

An  experimental  flow  bench  system  was  built  and  used  for  study  of  the  flow  de¬ 
cay  phenomena  under  Contract  AF04(61 1)-11620  (Ref.  2)  and  later  refurbished 
and  modified  for  use  under  Contract  F0461 1-68-C-0070  (Ref.  4).  For  the  cur¬ 
rent  program,  the  flow  bench  was  refurbished  again  and  a  number  of  changes 
were  made  to  the  system  to  satisfy  current  requirements.  Changes  included  in¬ 
stallation  of  larger  propellant  tanks,  an  agitator  in  the  main  propellant 
tank,  and  improved  instrumentation. 

A  schematic  diagram  of  the  flow  system  is  presented  in  Fig.  1  .  The  main  tank, 
or  run  tank,  and  the  catch  tanks  comprise  the  end  points  of  the  N2O4  flow 
path.  The  main  tank  is  a  10-gallon  flanged  spherical  stainless-steel  vessel 
contained  within  a  temperature-conditioned,  s tirred-water  hath .  The  catch 
tank  indicated  on  the  schematic  is  actually  two  5-gallon  stainless-steel 
cylinders  connected  in  parallel;  no  temperature  control  was  provided  for  these 
tanks.  Both  the  run  and  catch  tanks  are  connected  to  a  pressurization  and 
vent  system  which  supplies  gaseous  nitrogen  for  pressurization.  The  pressure 
control  system  maintains  both  the  tank  being  emptied  and  the  tanks  being 
filled  at  constant  pressure  during  any  flow  process.  The  controls  for  setting 
the  pressure  levels  are  operated  remotely. 

All  materials  of  construction  in  the  flow  bench  are  compatible  with  N^O^.  The 
tanks,  valves,  fittings,  and  fluid  lines  are  fabricated  of  type  304,  316,  321, 
or  347  stainless  steel.  Valve  seats  and  stem  packings  are  made  of  either 
Teflon  or  Kel-F.  Teflon  tape  is  used  as  the  thread  lubricant  on  pipe  fittings. 

Two  parallel  test  sections  are  provided.  The  flowrate  is  adjusted,  in  general, 
by  selecting  and  setting  the  desired  run  and  catch  tank  pressures.  However, 
each  flow  path  has  a  servo- operated,  remote-control  metering  valve  located 
downstream  of  the  test  section  which  can  be  used  for  additional  control.  For 
about  90  percent  of  the  tests,  this  valve  was  kept  full  open.  Between  each 
test  section  and  the  main  tank  is  a  heat  exchanger.  The  heat  exchangers  con¬ 
sist  of  a  length  of  1/4-inch  stainless-steel  tubing  in  an  open  water  bath. 

Each  heat  exchanger  bath  is  agitated  and  the  temperature  is  held  constant  by 
operating  an  on-off  valve  that  controls  the  flov,  of  coolant  water  from  a  re¬ 
frigeration  bath.  A  standpipe  returns  the  overflow  from  the  heat  exchanger 
bath  to  the  refrigeration  bath.  Ball  valves  are  provided  to  route  the  pro¬ 
pellant  flow  through  the  heat  exchangers  or  arouni  them  before  entering  the 
test  sections.  Also,  ball  valves  are  provided  to  select  either  the  run 
tank  standpipe  outlet  port  or  the  bottom  outlet  port,  depending  on  whether 
the  run  was  a  rapid  cooldown  test  or  a  long  or  slow  cooldown  test, 
respec  ti vely . 

After  completing  a  run,  propellant  can  be  discharged  to  a  waste  storage  tank 
or  recycled  to  the  main  tank  through  a  bypass  line  that  can  be  filtered  or 
not,  as  desired.  In  this  program,  the  recycled  propellant  was  never  fil¬ 
tered.  A  40-micron  filter  was  installed  in  the  system  feed  line  to  filter  all 
new  propellant  coming  into  the  flow  bench  from  the  storage  cylinders  (see 
F  1  C7  .  1  ) 
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Figure  1.  Schematic  Diagram  of  Experimental  Flow  Bench  System 


The  test  sections  consisted  of  small  stainless  filters  of  the  woven- 
wire  type.  These  filters  are  commercially  available  units  obtained  from 
Western  Filter  Company,  Inc.,  Los  Angeles  (basic  part  number  19310).  A  spe¬ 
cial  feature  of  this  filter  is  that  it  can  be  ?  '  *ed  to  or  removed  from  the 
system  with  minor  effort.  Thus,  it  was  parti<  ly  suited  to  the  needs  of 
this  program.  With  the  exception  of  one  filt  le  (a  SO-mesh  wire  screen), 
all  the  filters  used  in  this  program  were  in  tn-  ucron  range.  The  filter 
sizes  used  are  listed  below: 


Nil  cron 

Rating 

Size,  microns 

Nominal ^ 

Absolute  ^ 

2 

2 

10 

5 

5 

15 

10 

10 

25 

40 

40 

60 

50  mesh 

280 

(a) 

v  '95  to  98  percent  of  particles  above 
this  size  are  removed 

^100  percent  of  particles  above  this 
size  are  removed 

INSTRUMENTATION 

Nitrogen  tetroxide  flowrates  were  measured  in  each  of  the  parallel  test  sec¬ 
tions  by  Fischer-Porter,  radio- frequency -type  turbine  flowmeters.  Propel¬ 
lant  temperatures  were  measured  in  the  main  tank  and  in  the  line  between 
the  heat  exchanger  outlets  and  the  test  sections  (as  shown  in  Fig.  1  )  by 
iron-cons tantan  thermocouples  and  a  Pace  150-degree  reference  junction. 

Also,  thermocouples  were  installed  in  the  main  tank  bath  and  both  heat  ex¬ 
changer  baths  to  assist  in  setting  the  desired  test  conditions.  Pressures  in 
the  main  and  catch  tanks  were  determined  with  Taber  pressure  transducers. 
Pressure  drops  across  the  test  sections  were  measured  with  Data  Sensor  differ¬ 
ential  pressure  transducers.  These  flowrates,  temperatures,  and  pressures 
were  recorded  on  five  Hewlett-Packard  Mosely  dual -pen  strip  chart  recorders. 
Through  use  of  selector  switches,  heat  exchanger  bath  temperatures  could  be 
read  on  the  recorder  used  for  indicating  propellant  line  temperature,  but 
were  not  recorded  during  the  actual  test  runs. 

TEST  P  LAN 

A  test  plan  was  developed  and  approved  by  the  Air  Force  Project  Officer  be¬ 
fore  the  experimental  work  was  started.  It  was  designed  as  a  sequential 
plan  to: ( 1)  incorporate  features  to  minimize  the  required  number  of  tests, 

(2)  take  into  consideration  the  changeover  time  between  subsequent  tests, 
and  (3)  allow  flexibility  in  investigating  problem  areas  as  they  would  be 
identified  during  the  test  program. 

The  major  independent  variables  to  be  controlled  and  measured  were:  propel¬ 
lant  composition,  test  section  hardware,  initial  prope  1.'  ant  temperature, 
temperature  drop  before  reaching  the  test  section,  time  during  which  the 
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temperature  drop  is  imposed,  initial  flowrate,  and  saturation  condition 
(i.e.,  whether  the  1-204  was  used  as  received,  or  doped  with  a  small  amount 
of  iron  pentacarbonyl  to  ensure  iron  saturation) .  A  large  matrix  of  poten¬ 
tial  tests  was  prepared  with  a  preliminary  selection  of  the  combinations  of 
test  conditions  and  order  in  which  test  sets  would  be  conducted.  However,  it 
was  recognized  at  the  outset  that:  (1)  the  conditions  for  the  test  sets  to  be 
made,  (2)  the  number  of  tests  within  each  test  set  that  would  be  required  to 
obtain  the  necessary  information  (including  replications  to  help  resolve 
complex  and  unusual  behavior),  and  (3)  the  order  of  test  sets  would  differ 
somewhat  from  those  given  in  the  matrix  as  a  result  of  sequential  decisions 
(based  on  information  from  earlier  tests)  which  would  be  made  during  the 
actual  test  program.  This  preliminary  test  plan  outlined  many  more  poten¬ 
tial  tests  than  realistic  estimates  would  predict  could  be  performed  during 
the  program;  best  efforts  were  made  to  make  as  many  of  these  tests  as 
possible. 

During  the  program,  tests  were  made  with  variations  of  all  seven  of  the  major 
independent  variables  originally  selected.  These  variables,  and  the  levels 
or  ranges  of  values  employed  are  summarized  in  Table  I. 

TABLE  I.  INDEPENDENT  VARIABLES  IN  EXPERIMENTAL  TESTS 


1 

Variable 

Number  of  Levels 

Range  of  Values 

N204  Composition 

2 

Green  and  red-brown 

Test  Section  Hardware 

7 

2-micron  through  50-mesh 
filters 

Initial  Propellant 
Temperature,  F 

3 

50,  75,  and  100 

Temperature  Drop,  F 

7 

(approximate) 

0  through  34 

Time  for  Temperature 
Drop 

0 

Am 

Short  (during  fiow)  and  long 
(approximately  2  days) 

Initial  Flowrate,  gpm 

5 

(approximate) 

0.14  to  approximately  0.5 

Saturation  Condition 

2 

As  received,  or  doped  with  a 
small  amount  of  iron 
pentacarbonyl 

OPERATING  PROCEDURE 

The  run  tank  was  charged  with  nitrogen  tetroxide  from  the  appropriate  storage 
cylinder,  through  a  40-micron  filter  in  the  fill  line. 

For  certain  test  series,  a  contractual  requirement  was  to  use  propellant 
saturated  with  flow  decay  material.  To  ensure  that  saturated  conditions  ex¬ 
isted,  a  small  quantity  (an  amount  corresponding  to  "one  saturation  dose"  of 
iron  using  the  approximate  solubility  data  trom  Ket.  b)  of  iron  pentacarbonyl 
was  reacted  with  the  nitrogen  tetroxide  in  the  run  tank  to  foim  nitrosyl 
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tetranitratoferrate ,  the  flow  decay  compound.  The  iron  pentacarbonyl  was 
placed  in  a  microliter  syringe  with  a  24- inch  long  needle  and  injected 
directly  into  the  run  tank.  A  septum  and  ball  valve  arrangement  was  de¬ 
signed  and  mounted  on  the  top  flange  of  the  run  tank  for  this  purpose.  The 
tank  contents  were  agitated  during  this  addition. 

Both  long  and  rapid  cooldown  tests  were  conducted  in  this  program.  In  the 
"rapid"  cooldown  tests,  the  propellant  was  allowed  to  stand  in  the  run  tank 
for  a  minimum  of  2  hours  at  a  given  initial  temperature,  T  ,  before  a  test 
was  started.  The  main  tank  agitator  was  turned  on  periodically  to  help  equal¬ 
ize  the  temperature.  The  heat  exchanger  baths  were  then  adjusted  to  the  de¬ 
sired  temperature.  Final  prerun  preparations  consisted  of  setting  the  ball 
valves  in  the  flow  system  to  achieve  the  desired  flow  path.  If  the  propel¬ 
lant  was  to  be  subjected  to  a  temperature  drop  before  entering  the  test  sec¬ 
tion,  it  was  flowed  through  the  heat  exchanger  where  it  would  undergo  rapid 
cooldown.  If  the  propellant  was  bypassed  around  the  heat  exchanger,  i.e., 
no  temperature  drop  was  imposed,  the  run  would  be  a  2ero  temperature  drop 
test  but  still  a  part  of  the  rapid  cooldown  test  series. 

For  the  long  (or  slow)  cooldown  tests,  new  propellant  was  transferred  from 
the  storage  cylinder  to  the  flow  bench  and  allowed  to  stand  at  a  selected 
starting  temperature,  T  ,  for  at  least  2  hours.  Then  the  propellant  was 
cooled  slowly  over  a  period  of  about  48  hours  to  result  in  the  desired  tem¬ 
perature  drop.  The  propellant  was  maintained  at  the  final  temperature  for 
at  least  2  hours  to  ensure  that  temperature  equilibration  was  reached.  The 
run  was  then  conducted  with  the  propellant  undergoing  no  further  temperature 
changes  before  entering  the  test  section. 

To  begin  a  run,  the  main  tank  was  pressurized  to  the  selected  level  with  gas¬ 
eous  nitrogen.  The  catch  tank  was  then  pressurized  and  the  vent  regulator 
was  set  to  the  proper  pressure.  Flow  was  started  by  activatirg  the  remote- 
operated  ball  valve  at  the  catch  tank  inlet.  Flowrates  through  the  test  sec¬ 
tions  were  adjusted  as  necessary  by  making  changes  to  the  vent  regulator 
pressure  setting  and/or  the  main  tank  pressure  setting.  With  the  exception 
of  a  few  runs  (i.e.,  runs  with  very  low  flowrates)  the  metering  valves  were 
wide  open.  In  most  cases,  runs  or  a  series  of  runs  were  car'ied  to  pronel- 
lant  depletion.  At  this  point,  the  propellant  was  either  returned  to  the 
main  tank  through  the  bypass  line,  or  dumped  in  a  waste  storage  container. 

NITROGEN  TETROXIDE  COMPOSITIONS 

TVo  types  of  nitrogen  tetroxide,  red-brown  and  green,  were  used  in  this  pro¬ 
gram.  A  standard  nitrogen  tetroxide,  1-ton  shipping  container  was  emptied, 
cleaned,  and  sent  to  the  Western  Test  Range  (WTR)  and  another  to  the  Eastern 
Test  Range  (ETR)  to  be  filled  with  100  gallons  each  of  red-brown  N.,0^  and 
green  N'204,  respectively.  A  sampling  system  was  designed,  assembled,  and 
used  to  take  samples  for  chemical  analysis  from  each  of  the  two  propellant 
supply  tanks  after  they  were  returned.  These  samples  were  analyzed  by 
Rocketdyne  using  conventional  "wet"  techniques  (per  specification  MIL-P- 
26S39C,  and  in  some  cases  MI L-P- 26539B)  except  for  the  water  analyses,  which 
utilized  nuclear  magnetic  resonance  (*H  nmr)  techniques  with  N20,i-hl2n  stan¬ 
dards  carefully  prepared  and  sealed  in  precision  nmr  tubes  during  an  earlier 
Rocketdyne  program. 


6 


Information  also  was  obtained  from  ETR  and  WTR  on  their  chemical  analyses  of 
the  propellant  sources  from  which  the  propellants  shipped  to  Rocketdyne  were 
taken.  Information  from  KTR  showed  that  the  N2O4  in  the  ready  storage 
vehicle  (RSV)  used  for  the  2  June  1971  filling  of  Rocketdyne's  supply  tank 
was  loaded  on  20  and  21  August  1970.  The  WTR  chemical  analysis  results  were 
obtained  for  a  sample  taken  from  the  RSV  on  1  June  1971.  Results  of  chemi¬ 
cal  analysis  were  received  from  ETR  for  a  propellant  sample  taken  from  the 
RSV  used  to  fill  the  supply  tank  used  in  this  program.  The  ETR  sample  was 
taken  on  19  April  1971,  which  is  reasonably  close  to  the  time  (17  June  1971) 
that  the  supply  tank  was  filled.  The  results  of  these  chemical  analyses  are 
summarized  in  Table  II. 

TABLE  II.  CHEMICAL  ANALYSIS  OF  NITROGEN  TETROXIDE  USED  FOR  TESTS 


Component 

Red-Brown  N.,0^,  wt.  % 

Green  N-0.,  wt.  % 

2  4* 

NO 

H2° 

N'204  assay 

N0C1 

0.048(a')  ,  O.OSS^'1 

99. 9(b) 

O.S4(a),  0.63(C) 
0.060(a),  0.10(c) 
99 .4 Ca')  ,  99. O^-1 
0.03^  ,  0.02(c') 

^ Rocketdyne  analyses;  ^WTR  cr.r.i/ses;  *-c^ETR  analyses 

It  can  be  seen  that  there  is  fairly  good  agreement  between  results  of  the 
various  analyses. 

There  was  some  concern  about  the  implications  of  the  low  NO  concentration 
measured  for  the  green  N2O4  relati  e  to  the  lower  limit  of  the  current  MIL 
specification  (0.60  wt .  %) .  ETR  jotained  an  NO  content  of  0.63  wt.  %  com¬ 
pared  to  0.54  wt.  %  obtained  by  Rocketdyne.  The  two  analyses  were  made 
within  a  short  interval  of  time  and  therefore  both  can  be  considered  appli¬ 
cable  analyses.  An  average  of  the  two  points  gives  a  value  of  0.59  wt ,  % 
which  is  approximately  equal  to  the  lower  limit  of  the  current  MIL  specifi¬ 
cation.  All  of  the  other  component  concentrations  are  within  the  MIL-P- 
26S39C  limits. 


EXPERIMENTAL  RESULT'S  AND  ANALYSIS 


DERIVED  FLOW  DECAY  PARAMETERS 

The  ideal  experiment  to  measure  flow  decay  would  hold  the  driving  force  (pres¬ 
sure  differential)  across  the  test  section  exactly  constant.  Unfortunately, 
this  would  be  extremely  difficult.  In  this  program,  the  main  tank  and  catch 
tank  pressures  were  controlled  during  a  run  at  constant  levels  (as  discussed 
in  the  Flow  Decay  Test  section),  within  th6  operating  band  of  the  pressure 
regulators.  The  net  effect  on  pressure  drop  across  the  test  section  (H) ,  as 
the  resistance  across  the  test  section  increased  (i.e.,  flow  decay  was  occur¬ 
ring)  could  cause  H  to  decrease,  remain  constant,  or  increase.  TTiis  possible 
range  in  trend  is  a  result  of  the  complex  interaction  of  several  factors: 
the  actual  increase  in  resistance  across  the  test  section  caused  by  deposi¬ 
tion  of  flow  decay  material;  changes  in  the  pressure  drop  in  the  remainder  of 
the  flow  system,  as  the  flowrate  changes;  and  fluctuations  in  main  tank  and 
catch  tank  pressures,  within  the  normal  operating  bands. 

These  relationships  are  partly  illustrated  in  Fig.  2.  The  lowest  of  the 
three  curves  is  the  calibration  curve  (corresponding  to  a  completely  clean 
filter)  and  was  generally  known  from  the  calibration  measurements  at  the 
beginning  of  a  run.  Points  1  and  2  represent  values  of  Q  and  H  measured  at 
two  times  during  a  run.  In  this  case,  H2  is  shown  as  greater  than  Hj,  al¬ 
though  this  is  not  always  true  (as  discussed  in  the  previous  paragraph).  The 
two  dashed  curves  through  Points  1  and  2  represent  the  hydraulic  character¬ 
istic  curves  for  the  filter  with  the  degree  of  clogging  present  at  each  time 
but,  since  it  is  not  possible  to  "freeze"  the  filter  at  a  given  point,  it  is 
not  possible  to  determine  these  two  curves  experimentally. 

Several  methods  and  parameters  were  considered  to  reduce  the  basic  experimen¬ 
tal  data  to  a  common  basis,  in  order  to  permit  direct  comparisons  between 
runs  and  for  use  in  correlation  analyses.  The  final  choice,  and  most  mean¬ 
ingful  parameters  to  characterize  the  amount  and  rate  of  flow  decay  are  de¬ 
fined  as 


X  = 


(1) 


Y 


(2) 


where  t^  and  t2  are  times  corresponding  to  Points  1  and  2  during  a  run,  Qj 
is  the  measured  flowrate  at  time  tt,  and  Qo'  is  the  best  estimate  of  the 
value  of  Q2  if  H  were  held  constant  at  the  value  (i.e..  Point  a  in 
Fig.  2).  Therefore,  Y  is  the  best  estimate  of  the  rate  of  fractional  de¬ 
crease  in  flow  that  would  occur  if  the  driving  force  across  the  test  section 
remained  exactly  constant  during  a  test. 

Each  filter  test  section,  with  the  exception  of  those  used  on  Runs  1  throucrh 
39  and  160  through  165,  was  f’ow  calibrated  with  N2O4  to  be  used  in  that 
particular  test  series.  These  sets  of  calibration  data  points  were  fitted 
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very  well  by  the  equation 


H*  =  k(Q*)n  (3) 

where  the  value  of  n  for  various  filters  was  very  consistent  and  approximately 
equal  to  2.  The  value  of  k  did  vary  among  the  different  filters.  Typical 
pressure  drop  versus  flowrate  calibration  cur/es  for  several  different  filter 
test  sections  are  shown  in  Fig.  3.  In  Fig.  3,  to  distinguish  between  filters 
having  the  same  nominal  pore  size  but  different  filter  areas,  the  letter  "H" 
or  "F"  was  added  to  the  micron  rating. 


It  is  assumed  that  Eq.  3  also  represents  hydraulic  characteristic  curves  for 
partly  clogged  filters,  with  n  still  equal  to  2,  but  with  different  values 
of  k.  Equ  ‘ion  1  can  then  be  written  as 


X 


(4) 


This  equation  was  used  to  reduce  the  experimental  data  to  a  common  basis. 


Another  parameter  that  is  of  some  value  in  analyzing  flow  decay  data  was  de¬ 
fined  as 

w  , 

Q* 


(5) 


and  represents  the  flowrate  defect,  or  fractional  drop  in  flowrate  between  the 
new  filter  case  and  the  point  of  interest.  Values  of  W  at  various  stages  in 
a  filter's  history  can  be  compared,  e.g. 


W2  -  W1 
X2  "  tl 


It  can  be  shown  that  Y  and  Z  are  related  approximately  by  the  expression 

Q, 


-7  —  V 

i'"YQT* 


(6) 


(7) 


A  third  parameter  was  derived  for  use  in  examining  the  experimental  flow  decay 
data,  as  outlined  below.  The  flowrate  can  be  expressed  in  terms  of  a  dis¬ 
charge  coefficient,  C,  and  the  filter  area.  A,  as 


Q  = 


(8) 


The  discharge  coefficient  will  increase  as  the  filter  becomes  partly  clogged, 
and  is  a  function  of  many  variables  that  influence  the  deposition  of  flow 
decay  material,  including  a  filter  Reynolds  Number,  which  can  be  defined  as 


Rc 


DQP 

AU 


(9) 


If  it  is  assumed  that  C  is  proportional  to  the  free  filter  area,  Af,  (i.e., 
total  area  of  openings  available  for  flow)  at  any  time_.  and  that  the  amount 
of  effective  free  filter  area  that  is  clogged  by  flow  decay  material  (if  at 
all)  is  proportional  to  the  quantity  of  N2O4  per  unit  filter  area  that  has 
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flowed  through  the  filter,  then  the  following  equations  can  be  developed 


1 


C 


1 


C 

o 


Q(t)dt 


Q(t)dt 


L 


C 


2 


(10) 


(ID 


(12) 


where  Eq.  11  and  12  apply  for  constant  H  (i.e.,  C  is  the  same  for  both  t^ 
and  t2)  and  V^2  Is  tl'e  total  volume  of  N20,i  passing  through  the  filter  be¬ 
tween  time  tj  and  time  t2 .  Combining  Eq.  8  and  12,  plus  defining  0o  as  the 
new  filter  (i.e.,  calibration)  flowrate  results  in 


L 


(13) 


It  can  be  shown  that  X  and  L  are  related  approximately  by  the  expression 


2AX 

=  Q*  (t  2  -  tiy'(2  -  X) 


(14) 


ANALYSIS  OF  EXPERIMENTAL  ERROR 

Only  a  small  amount  of  flow  decay  occurs  with  the  type  of  filters  and  condi¬ 
tions  being  used  in  these  tests;  therefore,  a  number  of  experimental  error 
analyses  were  made  to  provide  guidance  in  the  data  reduction  work.  After  con¬ 
sidering  both  the  precision  of  various  experimental  measurements  and  the 
precision  associated  with  each  of  several  techniques  for  extracting  numbers 
from  the  recorder  charts,  it  was  possible  to  choose  a  convenient  technique 
giving  an  acceptable  precision  in  the  final  data  values  (i.e.,  flowrate  and 
pressure  drop  across  the  test  section  plus  corresponding  values  from  the 
filter  calibration  curves  at  various  times  during  a  run).  The  error  analy¬ 
sis  techniques  and  results  are  outlined  in  the  following  paragraphs. 

A  group  of  n  experiments  is  considered  in  which  pj  represents  the  true  value 
of  tne  outcome  of  the  i^  experiment.  It  is  recognized  that  there  is  experi¬ 
mental  error  at  each  step  and  thus  the  observed  values  correspond  to  a  ran¬ 
dom  variable  y^  such  that 


in  which  is  the  random  error  term.  It  is  assumed  that  the  are  all  in¬ 
dependent  and  have  expected  value  zero.  If  it  is  further  assumed  that  they 
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are  normally  distributed,  then  a  S5  percent  confidence  interval  on  p  is 

v.  i  1.96  o.  (IS) 

■  1  1 

where  a.  is  the  standard  deviation  of  v.  (and  hence  of  r.),  i.e., 

l  'i  i 

J.  2  =  var(v.  )  =  E(v.  -  E'vl)  2 
i  "l  'i 

=  E(e.2) 

Thus,  knowing  the  variance  of  y.  allows  establishment  of  error  bounds  on  the 
experimental  outcome. 

In  the  current  situation,  y^  is  a  function  of  several  independent  variables. 
Dropping  the  subscript  i  for  convenience,  this  can  be  written  as 

y  =  f(V-’Vm) 

where  v,,.. .,vm  are  the  values  of  the  independent  variables  (recorder  chart 
readings,  for  example),  and  which  are  subject  to  experimental  errors.  In 
general,  the  moments  of  y  (mean  and  variance)  are  functions  of  the  moments  of 
the  vj .  Computing  these  moments  can  be  a  very  arduous  task  unless  f  is  ex¬ 
tremely  simple.  In  many  cases,  this  can  not  be  done  in  closed  form;  there¬ 
fore,  it  may  be  necessary  to  resort  to  numerical  integration  or  possibly  the 
use  of  Monte-Carlo  techniques. 

However,  if  the  variances  of  the  Vj  are  small,  approximate  moments  of  y  can 
be  computed  by  expanding  f  in  a  generalized  Taylor's  series  about  the  means 
of  the  Vj .  Thus,  a  first-order  approximation  to  y  is  given  by 


y  =  f(u  ,...u  )  -  Z  §-  (x  -  u  ) 
J  =  1  )  J 


where  the  u-  are  the  means  of  v. 

J  J 


E(v.)  -  u. 


It  follows  that 


E(y)  =  f(Ul,  ...,U)  +  Z  E(v  -  u,) 


.  .  9v . 
1 


j  y 


var(y)  =  E(y  -  E(y)) 


t  ^  <v> '  ’J1] 


r  v  9f  nr 

£  ^  ITT”  37T"  "  ^nHv-  “  u;) 

£«1  j  =  l  3v£  9vj  x,  X,  j  ) 


If  it  is  assumed  that  v*  is  independent  from  v^,  then, 


E(v i  -  V(vj  -  >y  = 


JL  =  j 


0  8,  #  j 


Therefore,  the  final  result  is 


■  .i  m  v 

j=i  j  J 


(18) 


Comparing  Eq.  "3  with  Eq.  IS,  it  can  be  seen  that  the  problem  of  determining 
a  confidence  interval  for  the  experimental  outcome  has  been  reduced  to  esti¬ 
mating  the  variances  of  the  independent  variables. 

Equation  18  evaluated  with  X  as  the  dependent  variable  (i.e.,  y  =  X)  and 
using  the  expression  in  Eq.  4  to  define  f.  The  final  evaluation  of  this 
equation  was  made  with  the  following  values  for  the  means  and  variances  of 
the  experimental  variables. 


Independent 

Variable  (v.) 

J 

Nominal 

Mean  Value 

Es  timated 

2 

Variance  (a .  ) 

J 

j  Qj  and  Q2 

0.25  gpm 

2.5  x  10  7  gpm2 

Hj  and  H2 

7  psia 

-4  2 

4  x  10  psia 

C1  ‘  t2 

10  minutes 

-- 

tj  and  t2 

-- 

.  ,o-4  •  2 

1  x  10  min 

The  estimated  maximum  experimental  error  (95%  level)  in  each  case  is  two 
times  the  square  root  of  the  estimated  variance.  With  these  values,  and 
using  Eq.  18,  the  estimated  standard  deviation  of  X  is  0.0035  (i.e.,  a 
change  of  0.35%  in  flow).  This  estimated  value  is  essentially  independent 
of  the  time  interval  (t2  *  tj)  because  of  the  small  experimental  error  in 
measuring  time.  This  value  was  then  used  in  various  statistical  tests  ap¬ 
plied  to  the  data  to  determine  whether  flow  decay  occurred  or  not.  For  ex¬ 
ample,  experimental  runs  for  which  the  values  of  X  are  greater  than  1.282  0 
(i.e.,  0.0045)  can  be  said  to  represent  cases  of  flow  decay  (X  >  0)  at  the 
90  percent  confidence  level.  This  value  of  1.282  o  is  the  deviation  from 
the  mean  corresponding  to  a  confidence  level  of  90  percent  for  a  one-sided 
test  (normal  distribution) . 

CAT A  REDUCTION  AND  DERIVED  RESULTS 


Values  of  Q  and  H  were  taken  from  the  recorded  traces  for  each  run  at  several 
times:  (1)  as  near  to  the  beginning  of  the  run  as  practical,  (2)  at  the  end 
of  the  run.  and  (3)  at  one  or  more  intermediate  times,  especially  including 
any  inflection  points  within  a  run.  Generally,  the  slopes  of  the  Q  and  H 
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traces  were  fairly  constant  for  at  least  several  minutes  at  a  time.  A  com¬ 
puter  program  was  written  and  used  to  calculate  values  of  various  parameters 
(X,  Y,  Zj,  Zy,  W,  V,  L,  Re)  for  each  significant  time  interval  within  each 
run.  A  conplete  set  of  this  computer  output  is  given  in  the  Appendix.  These 
derived  results  were  then  used  with  a  variety  of  techniques  to  organize, 
generalize,  and  correlate  the  data. 

One  of  the  first  efforts  in  organizing  the  derived  results  was  to  develop  a 
large  matrix  showing  a  brief  indication  of  the  results  of  each  experimental 
run.  A  summary  of  results  for  all  the  flow  decay  tests  conducted  in  this  pro¬ 
gram  is  given  in  Tables  III,  IV,  and  V.  The  test  conditions  for  the  runs  are 
given  along  the  side  and  top  of  each  table  and  are  for  the  most  part  self- 
explanatory.  10  differentiate  between  filters  having  the  same  pore  size  but 
different  filter  areas,  the  letter  "H"  or  "F"  is  added  to  the  pore  size 
rating.  The  temperatures  listed  are  nominal  values;  the  actual  values  for 
each  run  are  given  in  Appendix  A. 

Each  box  in  these  tables  contains  information  about  the  tests  chat  were  run 
under  the  conditions  represented  by  the  values  of  the  seven  independent 
variables  corresponding  to  that  particular  box.  The  first  number  (or  group 
of  numbers)  in  each  box  is  the  test  run  number(s) .  The  second  entry  is  a 
simplified  indication  of  whether  or  not  that  particular  run  produced  flow 
decay.  The  word  "yes"  indicates  that  the  values  of  X  and  Y  were  high  enough 
to  indicate  that  flow  decay  occurred  (i.e.,  to  reject  the  hypothesis  that 
X  and  Y  are  equal  to  zero)  at  the  90  percent  level  of  confidence.  The 
error  analysis  and  estimated  standard  deviation  of  X  and  Y  as  outlined  in 
the  .Analysis  of  Experimental  Error  section  were  used  to  establish  the  accep¬ 
tance  region  for  the  above  test  at  the  90  percent  level  of  confidence.  The 
word  "no"  is  given  as  the  second  entry  in  a  l:ox  if  the  value  of  Y  is  low 

enough  to  indicate  that  no  significant  flew  decay  occurred  (for  this  purpose, 

the  minimum  level  of  significant  flow  decay  was  defined  to  be  Y  •  0.0$  %/min) 

at  the  90  percent  level  of  confidence.  The  abbreviation  "Unr"  is  given  if 

the  values  of  X  and  Y  do  not  permit  placing  the  result  in  one  of  the  first 
two  categories.  The  third  entry  in  each  box  is  the  value  (or  range  of 
values)  or  Y,  in  %/min.  In  some  instances,  where  tw©  or  more  runs  were  made 
under  a  given  set  of  test  conditions,  some  runs  did  and  other  runs  did  not 
exhibit  flow  decay.  In  these  instances,  the  table  box  is  divided  into  two 
or  three  sections  and  the  appropriate  information  is  given  for  each  group 
of  tests  in  the  format  described  above. 

There  are  many  observations  that  can  be  made  and  conclusions  that  can  be 
drawn  from  Tables  III,  IV,  and  V.  Many  of  these  are  incorporated  in  the  con¬ 
clusions  and  discussions  of  the  effects  of  the  independent  variables,  which 
are  given  in  the  remaining  sections  of  this  report.  However,  several  re¬ 
sults  should  be  noted  at  this  point,  including  the  following.  There  were 
no  cases  of  validated  flow  decay  with  filters  as  coarse  as  4U-micron. 

(This  was  also  true  for  50  mesh.)  There  were  five  runs  (numbers  1,  3,  4, 

10,  and  167)  in  which  the  results,  considered  without  the  other  experimental 
results,  would  indicate  valid  flow  decay.  However,  in  each  of  these  cases, 
there  are  other,  more  compelling  reasons  to  conclude  that  true  flow  decay 
will  not  occur  under  those  conditions.  For  each  of  these  five  anomalous 
runs,  there  are  much  larger  numbers  of  data  points  at  the  same  conditions 
(and  others  with,  e.g.,  higher  temperature  drops  and,  smaller  filter  pore 
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TABLE  III.  SUMMARY  OF  FLOW  DECAY  TESTS  USING  IRON-SATURATED  N„0 
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ni ti al  Propellant  Temperature 
Nominal  Temperature  Drop, 


TABLE  V.  SUMMARY  OF  FLOW  DECAY  TEST  RESULTS  USING  GREEN  AS-RECEIVED  N,0 


sizes)  for  which  flow  decay  was  not  present.  Therefore,  it  is  judged  valid 
to  conclude  that  flow  decay  will  not  occur  under  any  of  the  conditions  tested 
with  filters  of  nominal  pore  size  at  least  as  large  as  40  microns. 

There  are  other  trends  that  can  be  observed  merely  from  examination  of  these 
tables.  The  incidence  of  flow  decay  (and  the  rate  of  decay)  is  increased 
by  increasing  the  temperature  drop  and  by  decreasing  the  filter  pore  size. 
Although  some  additional  trends  may  be  suspected  from  examination  of  these 
tables,  caution  should  be  exercised  since  the  situation  is  more  complex  than 
it  may  appear  at  first.  There  are  other  independent  variables  that  are  not 
explicitly  listed  in  Tables  III  through  V.  These  include  initial  flowrate 
(which  was  measured,  but  is  not  shown  in  these  tables)  and  the  "history 
parameter,"  V,  which  expresses  the  total  amount  of  propellant  per  unit  filter 
area  which  has  passed  through  the  filter.  Another  complexity  arises  because 
there  appear  to  be  many  interactions  among  the  various  independent  variables. 

Part  of  the  project  was  devoted  to  going  beyond  the  yes-no  level  of  informa¬ 
tion  about  the  incidence  of  flow  decay,  and  to  perform  correlations  to  es¬ 
tablish  at  least  a  preliminary  basis  for  representing  the  amount  of  flow  de¬ 
cay  that  occurred  under  various  experimental  conditions.  Because  of  the 
large  number  of  independent  variables  and  complexity  of  the  flow  decay  phen¬ 
omenon,  it  was  necessary  to  rely  on  other  statistical  plus  multiple  regres¬ 
sion  analyses.  These  are  described  in  the  next  section. 

MULTIPLE  REGRESSION  ANALYSIS 

Multiple  regression  techniques  were  used  to  analyze  the  data.  This  type  of 

analysis  provides  -  'onvenient  tool  to  determine  if  various  independent 

variables  have  a  it  effect  on  the  outcome  of  a  test.  This  tech¬ 
nique  is  also  va  Jetecting  cross  effects  or  interactions  between 

the  independent  '  In  particular,  a  step-wise  multiple  regression 

was  performed,  in  arious  combinations  of  the  independent  variables 

were  allowed  in  the  t  .  O.  ;ssion,  and  selections  were  made  to  delete  variables 
that  were  of  low  value  in  tne  regression  and  to  add  (or  retain)  variables 
that  were  of  high  value  (i.e.,  contributed  significantly  to  reducing  the  re¬ 
sidual  sum  cf  squares). 


To  provide  some  guidance  for  the  regression  analyses,  theoretical  analyses 
were  made  of  the  physical  processes  that  may  be  occurring  during  flow  decay. 
However,  since  the  actual  complex  mechanisms  of  flow  decay  are  only  partly 
understood,  and  are  generally  not  very  amenable  to  analysis,  it  was  necessary 
to  perform  most  of  the  regression  efforts  in  an  empirical  manner.  Approxi¬ 
mately  100  sets  of  step-wise  multiple  regression  analyses  were  made;  each 
case  comprised  a  series  of  approximately  5  to  50  steps  in  which  various  com¬ 
binations  of  the  independent  variables  were  permitted. 


The  nine  primary  independent  variables  considered  in  the  regressions  were: 
Reynolds  number  (Re);  nominal  filter  pore  size  (D) ;  temperature  drop  (AT); 


initial  ten^ierature  (T  )  ;  the  average  flowrate  per  unit  filter  area  (Q/A)  ; 
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per  unit  filter  area  that  passed  through  the  filter  during  the  run;  and 


three  2- level  parameters.  The  2-level  parameters  were:  (1)  the  type  of 
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N2O4  (green  or  red-brown),  (2)  whether  or  not  the  propellant  was  doped  with 
a  small  amount  of  iron  pentacarbonyl  to  ensure  saturation  or  was  tested  as- 
received,  and  (3)  whether  the  tenperature  drop  was  imposed  rapidly  (during 
flow)  or  slowly  (over  approximately  2  days).  The  two  possible  levels  for 
each  of  these  three  parameters  were  assigned  values  +1  and  -1,  respectively, 
for  use  in  the  regression  analyses.  In  addition  to  these  nine  primary  in¬ 
dependent  variables,  a  large  number  of  combinations  of  these  variables  were 
also  allowed  as  possible  independent  variables. 

Regression  analyses  were  made  for  the  total  set  of  test  data  and  for  many 
partitioned  subsets  of  the  data.  There  are  over  1,000,000  combinations  of 
the  nine  primary  variables  with  the  nominal  number  of  levels  tested  for  each 
during  this  program.  Therefore,  it  can  be  readily  seen  that  it  was  impor¬ 
tant  to  examine  the  data,  where  possible,  for  subsets  in  which  there  are 
fewer  independent  variables  and  hence,  fewer  possible  combinations  of  condi¬ 
tions.  The  subsets  included  such  groupings  as;  all  runs  for  which  there 
was  significant  flow  decay;  all  runs  with  green  N2O4  for  which  there  was 
significant  flow  decay;  all  runs  with  green  N2O4,  tested  as-received,  sub¬ 
jected  to  a  rapid  temperature  drop;  etc.  The  multiple  regression  results  for 
the  subsets  were  more  useful  than  those  for  the  total  set  of  data  in  deducing 
effects  and  interactions  of  the  independent  variables. 

To  demonstrate  the  ability  of  multiple  regression  analysis  techniques  to 
extract  information  from  complex  data,  the  following  simple  example  from  the 
multiple  regression  analyses  performed  during  this  project  is  discussed  in 
some  detail.  Several  analyses  were  performed  for  a  partitioned  subset  of 
the  test  data  consisting  of  Runs  120  through  128  and  130  through  134.  For 
these  runs,  the  filter  pore  size  and  initial  temperature  were  held  constant 
at  5u  and  50  F,  respectively.  Moreover,  the  propellant  for  all  runs  was 
red-brown  N2O4,  saturated  with  respect  to  iron,  and  was  cooled  rapidly.  The 
only  parameters  that  were  allowed  to  vary  were  AT,  V,  Q/A,  and  T  .  It  had 
been  found  from  previous  analyses  that  inclusion  of  Re  did  not  improve  the 
regression  as  compared  with  the  use  of  Q/A  and  D  alone.  For  these  particular 
runs,  V  and  Q/A  were  highly  correlated  and  therefore  V  was  also  dropped  as 
an  independent  variable.  It  then  follows  that  the  rate  of  flow  decay,  Y, 
is  some  function  of  only  AT  and  Q/A  plus  experimental  error.  Typically  in  a 
regression  analysis,  it  is  assumed  that  the  dependent  variable  (or  some 
transformation  of  it)  is  normally  distributed  with  constant  variance  but 
with  mean  depending  on  the  independent  variables,  i.e., 

E (Y)  =  f (AT,  Q/A)  var(Y)  =  a2  (19) 

where  E  =  expectation  and  var  =  variance.  Assuming  that  f  is  reasonably 
well  behavod,  it  can  be  approximated  by  an  expansion  about  nominal  values 
of  AT  and  Q/A.  In  this  case,  the  expansion  to  be  considered  is 

f (AT,  Q/A)  =  ax  +  a2  AT  ♦  a3(Q/A  -  1.5)  +  &4  AT (Q/A  -  1.5)  (20) 

The  a^  are  the  unkncwn  regression  coefficients  and  are  to  be  estimated  using 
the  standard  least  squares  criterion.  Under  the  model  assumptions  given, 
these  estimates  will  be  optimum  ir.  the  sense  that  they  have  minimum  mean 
squared  error  loss  among  all  possible  invariant  estimates  (those  estimates 
invariant  under  change  of  scale  or  location  of  origin).  Th'iy  are  also  mini¬ 
mum  variance  estimates  among  all  unbiased  ones. 
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Equation  20  car  be  used  for  several  purposes.  First  of  all,  the  hypothesis 

H  :  a.  =  0 
o  1 

can  easily  be  tested  against  the  alternatives 

H,:  a.  4  0,  H,:  a.  >  0,  HT:  a.  <  0 

1  i  2  i  3  i 

using  the  standard  F-statistic.  The  confidence  at  which  this  hypothesis  can 
be  rejected  gives  a  measure  of  the  confidence  with  which  it  can  be  concluded 
that  the  parameter  for  which  is  a  coefficient  has  an  effect  upon  flow 
decay.  The  sign  of  the  estimate  of  a.  gives  the  indication  of  the  direction 
of  this  effect. 

From  the  regression  analysis,  the  estimates  of  the  coefficients  aj ,  a2,  a3, 
a4  are  0.12,  0.025,  0.093,  and  -0.003  respectively.  The  coefficient  a2  is 
considered  to  represent  the  main  effect  of  the  variable  AT  (defined  here  as 
the  slope  of  the  regression  curve  at  the  nominal  values  of  AT  and  Q/A) . 

Since  a2  >  0,  the  regression  equation  predicts  that  a  temperature  drop  will 
tend  to  enhance  flow  decay.  The  standard  deviation  of  this  estimate  is 
given  as  0.013.  Thus,  it  can  be  said  that  this  effect  is  real,  with  a  high 
degree  of  confidence;  in  fact,  the  hypothesis  that  a2  =  0  can  be  rejected 
at  a  92  percent  confidence  level.  Similarly,  the  value  of  a3  as  0.093  in¬ 
dicates  that  (for  AT  =  0)  the  effect  of  Q/A  is  to  enhance  flow  decay.  The 
confidence  with  which  a3  =  0  can  be  rejected  is  again  quite  high  (93  percent) . 

In  this  example,  a  cross  effect  between  AT  and  Q/A  was  allowed  through  in¬ 
clusion  of  the  term  a4  AT(Q/A  -  1.5).  The  value  of  -0.003  for  a4  indicates 
that  the  effect  of  AT  is  decreased  as  Q/A  is  increased.  This  follows  from 
the  fact  that 
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3 AT  a2  +  a4 

=  0,025  -  0.003  (Q/A  -  1.5)  (21) 

Similarly,  the  effect  of  Q/A  is  decreased  as  AT  is  increased.  The  standard 
deviation  of  a4  is  0.008.  The  hypothesis  that  a4  =  0  can  be  rejected  only 
at  a  28  percent  level  of  confidence.  This  is  quite  low;  therefore,  caution 
should  be  exercised  in  drawing  any  general  conclusions  concerning  the  ex¬ 
istence  of  the  negative  cross  effect  between  AT  and  Q/A. 

A  plot  of  the  data  points  and  regression  curve  predictions  for  this  case 
is  given  in  Fig.  4.  The  positive  Q/A  effect  for  any  given  AT  is  seen  by 
noting  that  the  family  of  curves  all  have  a  positive  slope.  The  positive 
AT  effect  is  seen  by  noting  that  each  curve  in  the  family  has  a  higher  AT 
indexing  it  than  the  one  immediately  below  it.  Finally,  the  small  cross  ef¬ 
fect  between  AT  and  the  Q/a  variable  is  evidenced  by  the  changing  slope  of 
the  three  curves  plotted. 

Results  of  several  of  the  more  useful  multiple  regression  analysis  cases  for 
partitioned  subsets  of  the  data  are  given  in  Tables  VI  and  VII.  The  case 
numbers  giver,  in  these  tables  do  not  mean  that  these  were  the  first  eight 
analysis  cases  that  were  run. 
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Figure  4.  Regression  Equation  Predictions  and  Data  Points  (Regression 
Case  2  in  Table  VI) 


TABLE  VII.  FINAL  REGRESSION  EQUATIONS  FOR  SELECTED  DATA  SUBSETS 
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For  a  given  case,  a  single  entry  in  a  column  of  Table  VI  indicates  tha*  the 
variable  in  that  column  was  held  fixed  at  the  value  shown  for  all  tests  in¬ 
cluded  in  that  subset.  (Conversely,  all  test  data  for  the  specified  condi¬ 
tions  were  included  in  the  subset.)  If  three  numbers  are  given,  this  indi¬ 
cates  that  the  variable  was  included  in  the  set  of  independent  variables.  The 
top  entry  is  the  estimate  of  the  regression  coefficient,  the  second  is  the 
estimate  of  its  standard  deviation,  and  the  third  is  the  level  of  confidence 
at  which  the  hypothesis  can  be  rejected  that  the  regression  coefficient  is 
zero.  For  example,  in  Case  S,  the  estimate  of  the  regression  coefficient  for 
the  variable  GB  is  0.41,  its  standard  deviation  is  0.07,  and  there  is  99.9 
percent  confidence  that  this  coefficient  is  not  zero. 


The  value  of  a  given  coefficient  cart  be  interpreted  as  the  effect  of  the 
independent  variable  which  it  represents.  More  precisely,  the  effect  is 
the  slope  of  the  regression  curve  evaluated  at  the  nominal  values  0=7 
microns,  AT  =  15  F,  TQ  =  75  F,  GB  =  Q,  SAT  =  0,  RL  =  0,  V  =  15  gal/sq  in., 
and  Q/A  =  1.5  gal/min-sq  in.  (except  for  Case  2  for  which  the  nominal  AT 
value  used  was  zero).  In  some  instances,  there  is  an  extra  term  added  to 
the  coefficient.  For  example,  the  coefficient  for  AT,  Case  5,  is  given  as 
0.020  -0.008  GB.  This  is  to  be  interpreted  as  (0.020  -  0.O08)  =  0.012 
for  green  N2O4  and  (0.020  +  0.008)  =  0.028  for  brown  N2O4.  me  estimate  of 
the  standard  deviation  and  confidence  level  apply  to  the  0.020  portion 
only.  Evaluation  of  the  standard  deviation  corresponding  to  the  values 
0.012  and  0.028  would  involve  an  expression  requiring  the  variances  and 
covariances  of  the  AT  and  AT  •  GB  terms,  that  is 


a 

a 


2 


2 


2cov(a1,  a3) 


(22) 


where  aj  is  the  coefficient  of  AT,  a2  is  the  coefficient  of  AT  for  GB  =  1 , 
and  a,  is  the  coefficient  of  AT  •  GB .  These  computations  were  not  per¬ 
formed,  in  general. 


The  conclusions  that  were  drawn  from  all  of  the  regression  analyses  are  dis¬ 
cussed  in  the  final  subsection.  Effects  of  Independent  Variables,  which  fol¬ 
lows  the  considerations  of  threshold  effects. 


THRESHOLD  EFFECTS 

It  has  been  observed  that  there  often  appears  to  be  a  "threshold"  effect 
for  flow  decay  i.e.,  an  identifiable  boundary  between  a  range  of  the  varia¬ 
bles  for  which  no  flow  decay  occurs  and  a  range  for  which  flow  decay  occurs 
in  varying  amounts.  Such  thresholds  are  probably  not  sharp,  and  certainly 
depend  on  the  interactions  of  the  independent  variables.  Investigations 
of  threshold  effects  were  made  during  this  program.  Descriptions  of  gen¬ 
eral  threshold  behavior,  techniques  for  deducing  threshold  values  from  the 
experimental  data,  and  examples  of  threshold  effects  are  described  in  this 
subsection. 

The  experimental  model  assumed  in  Eq .  19  s  capable  of  describing  a  broad 
range  of  real  situations,  its  application,  however,  depends  upon  the 
ability  to  characterize  the  functional  form  of  E(Y)  in  terms  of  parameters 
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(regression  coefficients),  which  can  be  estimated  from  the  data.  In 
the  present  situation,  however,  very  little  is  known  from  physical  grounds 
concerning  the  behavior  of  flow  decay  as  a  function  of  the  various  inde¬ 
pendent  variables  that  have  been  defined.  In  the  example  presented  in 
the  Multiple  Regression  Analyses  subsection,  it  was  assumed  that  f  is  an 
analytic  function  of  its  variables  and  hence  can  be  approximated  by  a  trun¬ 
cation  of  a  generalized  Taylor's  series  expansion  about  the  nominal  values 
of  those  independent  variables.  This  expansion  was  assumed  to  be  valid 
only  in  a  region  where  there  was  positive  flow  decay.  The  reason  for  this 
restriction  is  due  to  the  possible  threshold  effect  for  flow  decay.  If 
Y  is  used  as  the  dependent  variable  denoting  flow  decay,  v  =  (v^,...,vn) 
denotes  the  vector  of  test  conditions  and  c  denotes  the  region  of  n-space 
for  which  there  is  positive  flow  decay,  then  the  threshold  that  is  being 
discussed  is  the  boundary  of  c.  In  particular,  the  expectation  of  Y, 
called  f  in  Eq.  19,  can  be  written  as 

(0  v  e  c' 

f(v)  =  _  (23) 

I  g(v)  V  E  C 

in  which  g  is  a  positive  function  and  c'  is  the  complement  of  c.  The  test 
data  appear  to  bear  out  the  existence  of  the  region  c'  (region  of  no  flow 
decay)  and  hence  Eq .  23  is  probably  an  accurate  description.  Clearly,  f 
cannot  be  approximated  by  a  series  expansion  for  all  v  because  it  is  not 
analytic.  However,  it  is  assumed  to  be  continuous.  This  continuity  assump¬ 
tion  makes  it  possible  to  estimate  the  boundary  of  c  from  estimates  of  g 
by  finding  those  v  such  that  g(v)  =  0,  thus  the  boundary  of  c,  3c,  is  given 
by 


3c  =  (v|g(v)  =  0}  (24) 

It  is  actually  the  function  g  that  is  being  approximated  by  a  truncated 
series  expansion.  The  coefficients  of  this  expansion  are  estimated  in  the 
regression  analysis  previously  discussed.  Since  g  is  only  defined  on  c, 
those  points  not  in  c  must  be  deleted  from  the  regression  analysis.  The 
problem  here,  however,  is  that  the  boundary  of  c  is  not  known  a  priori; 
therefore,  it  is  difficult  to  decide  on  the  basis  of  the  value  of  a  partic¬ 
ular  v  whether  or  not  that  data  point  should  bo  deleted  from  the  regression 
analysis.  The  problem  is  circumvented  by  using  only  those  points  for  which 
there  is  probable  true  flow  decay  (i.e.,  Y  >  0  at  a  specified  ievel  of  con¬ 
fidence)  and  hence  it  is  probable  that  v  e  c  (i.e.,  v  is  a  member  of  the 
set  c). 

Two  examples  of  this  threshold  phenomenon  are  described.  The  first  is  Case 
5  in  Table  VI.  For  this  case,  the  initial  propellant  tempera ture  was  held 
fixed  at  7S  F;  the  N-jO^  was  saturated  and  cooled  rapidly.  These  restric¬ 
tions  leave  14  of  the  original  26  data  points  in  this  case.  Of  these,  2 
had  a  filter  diameter  of  10  microns;  the  remainder  had  a  5-micron  filter. 

It  therefore  seems  useful  to  analyze  the  case  for  D  =  5  microns  because 
this  is  where  the  highest  confidence  in  the  estimated  regression  curve 
can  be  placed.  The  following  variables  are  left  to  vary:  temperature 
drop  (at),  history  parameter  (V),  and  average  flow  per  filter  area  (Q/A) . 
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The  regression  equation  for  Case  5  thus  reduces  to 

Y  =  9.035  +  0.0228  (AT  -  IS)  +  0.0008  (AT  -  IS) 2  -  0.02S0  (V  -  IS) 

♦  1.526  (Q/A  -  l.S)  -  0.0010  (AT  -  IS)  (V  -  15) 

♦  0.0176  (AT  -  15)  (Q/A  -  l.S)2  (25) 

where  D  is  set  equal  to  S  microns  and  GB  is  set  equal  to  1  (i.e.,  green 
N1O4)  .  Those  values  of  AT,  V,  and  Q/A  for  which  Y  in  Eq.  25  is  zero  will 
thus  approximate  a  portion  of  the  boundary  of  c  when  the  variables  take  on 
the  previously  prescribed  values. 

In  practice,  however,  this  approximation  is  only  valid  when  Eq.  25  repre¬ 
sents  a  reasonably  good  fit  of  the  data.  Further,  since  the  values  of  AT, 

V,  and  Q/A  for  which  Y  is  zero  will  lie  outside  the  region  for  which  data 
was  used  to  estimate  the  coefficients  of  Eq.  25,  then  extrapolation  errors 
may  be  introduced.  These  will  be  minimized  for  those  cases  in  which  the 
data  used  are  "close"  to  the  threshold.  This  appears  to  be  the  case  for 
the  present  example. 

Figure  5  is  a  plot  of  those  data  points  used  to  estimate  the  regression 
curve  for  which  the  N2O4  was  green  and  the  filter  diameter  was  5  microns. 
Further,  only  those  points  for  which  AT  was  less  than  or  equal  to  20  F 
are  displayed.  Also  shown  are  representatives  of  the  regression  equations. 

To  avoid  constructing  a  four -dimensional  graph  (Y  *  f(AT,  V,  Q/A)),  the 
following  mechanism  was  used  to  plot  Y: 

1.  The  vertical  axis  is  the  rate  of  flow  decay,  Y 

2.  The  horizontal  axis  is  the  average  flow  per  filter  area,  Q/A 

3.  The  family  of  curves  is  indexed  by  the  temperature  drop,  AT 

4.  The  history  parameter  was  mapped  onto  the  Q/A  coordinate  through 
linear  or  quadratic  mappings 

Step  4  may  require  further  comment.  The  line  for  AT  -  0,  e.g.,  is  the  re¬ 
gression  line  which  would  ideally  pass  through  both  AT  =  0  data  points 
(triangles)  in  the  case  of  a  perfect  curve  fit.  The  lower  AT  *  0  point  has 
associated  with  it  a  V  of  10.0  and  a  Q/A  of  0.88.  The  upper  point  has  a 
value  of  V  =  26. S  and  a  value  of  Q/A  =  2.01.  The  straight  line  in  the  V, 

Q/A  plane  given  by 

V  =  10.0  ♦  14.6  (Q/A  -  0.88) 

passes  through  these  two  points.  The  straight  line  labeled  AT  =  0  in 
Fig.  5  is  then  obtained  by  substituting  Eq.  26  into  Eq.  25  as  well  as 
setting  AT  =  0  in  Eq.  25,  The  result  is  a  linear  relationship  between  Y  and 
Q/A.  Equivalently,  the  AT  =  0  curve  in  Fig.  5  may  be  thought  of  as  a  pro¬ 
jection  onto  the  Y,  Q/A  plane  of  a  line  on  the  regression  surface  in  the 
three-dime  ns i oral  space  ,(Y,  V,  Q/A)  (see  Fig.  6).  This  line  on  the  regres¬ 
sion  surface  (line  AC  in  Fig.  6)  connects  the  two  points  which  are  given  by 

th»  on  Bstimatps  of  fhf*  1ow*»r  and  imn«r  AT  =  0  no-infs  ("no-infs  a 

0  ------  -  *  -  -  —  -  -  •  -  -  -  ~  c  ■»«•■■■■- 

and  b  in  Fig.  6).  In  general,  there  are  an  infinite  number  of  lines  con¬ 
necting  these  two  points  ard  which  lie  on  the  regression  surface.  The  one 
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chosen  here  is  that  which  projects  as  a  straight  line  onto  the  V,  Q/A  plane 
( line  AB  in  Fig .  6) . 

The  line  plotted  in  the  V,  Q/A  plane  (line  DE  in  Fig.  6)  can  \  '  seen  to  in¬ 
tercept  Y  =  0  for  Q/A  near  0.75.  This  intercept  is  the  threshold  value  of 
Q/A  predicted  by  the  regression  surface  for  V  =  7.7.  It  can  thus  be  seen 
that  the  plotting  procedure  used  for  Fig.  5  will  depict  thresholds  for  Q/A. 

The  information  not  represented  in  Fig.  5  is  the  value  of  V  at  which  these 
Q/A  thresholds  lie. 

The  AT  =  5  curve  has  three  data  points  associated  with  it.  These  three  points 
determine  a  quadratic  relationship  between  V  and  Q/A.  The  AT  =  10  curve 
again  uses  a  linear  transformation.  For  AT  =  20,  there  are  three  points  and 
hence  a  quadratic  function  would  suffice.  However,  it  was  easier  to  picture 
the  regression  curves  by  using  two  separate  linear  mappings.  The  (a)  curve 
is  meant  to  approximate  the  AT  -  20  points  (diamonds)  at  Y  =  0.36  and  Y  =  1.35, 
while  the  (b)  curve  approximates  those  at  Y  =  0.10  and  Y  =  1.35. 

Admittedly,  this  procedure  of  graphing  the  data  and  regression  curves  in  two- 
dimensions  is  somewhat  artificial.  However,  it  does  serve  two  extremely 
important  functions.  First,  it  shows  that  the  data  are  reasonably  well  fit 
by  the  multip’e  regression  equation  (Eq.  25).  Secondly,  it  provides  a 
visual  check  on  how  well  an  extrapolation  of  the  regression  curve  to  zero 
might  truly  represent  the  thresholds  of  the  independent  variables. 

Figure  7  represents  the  functional  relationships  between  AT,  V,  and  Q/A  at 
the  threshold.  These  are  found  by  setting  Y  =»  0  in  Eq.  25.  It  can  be  seen 
that,  for  each  given  V,  Q/'A  is  nearly  a  monotonical lv  decreasing  function  of 
AT.  Conversely,  for  each  V,  AT  is  a  decreasing  function  of  Q/A.  That  is, 
as  the  flow  per  filter  area  is  increased,  it  is  necessary  to  go  to  to  smaller 
temperature  drops  to  avoid  flow  decay.  The  lack  of  pure  monotonicity  may  be 
due  to  a  modelling  error  in  not  defining  the  exact  types  of  interactions 
possible.  Nevertheless,  the  goodness  of  fit  depicted  in  Fig.  5  indicates 
that  the  thresholds  presented  in  Fig.  7  are  reasonable  approximations. 

Case  1  in  Table  VI  is  considered  for  a  second  example  of  the  threshold  phen¬ 
omenon  for  N2O4  flow  decay,  as  predicted  by  a  regression  curve.  Attention 
is  restricted  to  that  subset  of  the  data  used  in  Case  1  that  employed  a 
rapid  cooldown  for  the  temperature  drop.  Figure  8  is  a  graph  of  these  data 
as  well  as  the  regression  surface  given  in  Table  VII,  and  evaluated  for 
rapid  cooldown. 

Plane  ABCD  is  a  three-dimensional  representation  of  this  regression  equa¬ 
tion,  giving  flow  decay  Y,  as  a  function  of  V  and  Q/A.  In  the  group  of 
runs  used  for  this  case,  the  following  variables  did  not  vary:  AT  =  20  F, 

T0  =  75  F,  D  =  5  micron,  and  use  of  as-received  (ASIS)  brown  N2O4.  Line 
BC  lies  in  the  plane  of  the  V  and  Q/A  axes  and  represents  the  intersection 
of  the  regression  plane  with  the  Y  =  0  plane.  Similarly,  lines  AD  and  CD 
lie  in  the  planes  formed  by  the  V  and  Y  axes  and  the  Q/A  and  Y  axes,  re¬ 
spectively,  and  represent  the  intersections  of  the  regression  plane  with 
these  plar.es.  Although  the  plane  formed  by  the  regression  curve  is  shown 
in  Fig.  8  as  ABCD,  the  plane  extends  beyond  AD. 
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Figure  S.  Example  of  Regression  Plane  and  Threshold  Effects  (Regression  Case  i) 


The  regression  equation  predicts  that  no  flow  decay  will  occur  if  the  data 
point  has  V,  Q/A  coordinates  that  fall  between  )ine  BC  and  the  origin  of  the 
axes.  These  values  of  V,  Q/A  are  below  the  threshold  where  the  regression 
curve  predicts  flow  decay.  All  other  positive  values  of  V  and  Q/A  are  at 
or  above  the  threshold  values  and  therefore  flow  decay  would  be  expected.  1 

EFFECTS  OF  INDEPENDENT  VARIABLES 

There  are  a  number  of  trends  and  effects  of  the  independent  variables  (all 
of  which  represent  or  are  derived  directly  from  the  test  conditions)  which 
can  be  deduced  from  the  test  data  and  the  results  of  the  regression  analy¬ 
ses  performed  on  the  data.  The  major  effects  are  discussed  in  this  sub¬ 
section.  The  reader  should  continually  keep  in  mind  that  extreme  caution 
must  br  taken  in  extrapolating  any  effects  and  trends  to  different  flow 
situations  and  conditions.  The  flow  decay  phenomenon  is  very  complex  and 
is  affected  by  a  large  number  of  independent  variables.  The  complexity 
is  increased  by  the  substantial  interaction  effects  between  the  variables. 

It  is  readily  apparent  that  the  effect  of  the  parameter  AT  is,  in  general, 
positive.  That  is,  flow  decay  is  increased  as  the  temperature  drop  as  in¬ 
creased.  This  conclusion  follows,  e.g.,  from  regression  Cases  2,  5,  and  7 
(Table  VT'  for  which  the  estimate  of  the  AT  coefficient  is  more  than  two 
standard  deviations  from  zero.  Case  4  also  lends  credence  to  this  conjec¬ 
ture  of  a  positive  AT  effect  although  the  confidence  level  is  only  0.62. 

The  negative  coefficients  for  Cases  3  and  8  are  not  significant  because  they 
fall  well  within  the  range  of  estimated  errors.  Case  6  is  a  little  more 
subtle.  The  data  indicate  that  the  slope  is  positive  (0.002  *  0.011  =  0.013) 
for  saturated  conditions,  while  negative  (0.002  -  0.011  =  -0\009)  for  un¬ 
saturated  conditions.  If  the  standard  deviation  of  the  -0.009  estimate  is 
computed,  however,  it  is  found  to  be  0.01.  Thus,  again,  the  negative  coef¬ 
ficient  could  be  attributed  to  data  scatter.  Conversely,  the  positive  coef¬ 
ficients  of  Cases  2,  5,  and  7  cannot  be  attributed  to  data  scatter  and 
hence  the  conjecture  of  a  positive  AT  effect  is  upheld. 

The  filter  pore  size  effect  on  the  other  hand  is  generally  negative,  mean¬ 
ing  that  flow  decay  is  decreased  as  the  pore  size  of  the  filter  is  in¬ 
creased.  This  conclusion  follows  from  Cases  5,  7,  and  8  for  which  the 
hypothesis  that  the  diameter  coefficient  is  zero  can  be  rejected  with  0.98, 
0.83,  and  0.28  levels  of  confidence,  respectively.  Admittedly,  the  0.28 
level  is  very  low  but  the  trend  is  still  apparent. 

A  possible  anomaly  in  the  effect  oi  D  may  appear  from  close  examination  of 
Case  5.  The  estimate  of  the  diameter  effect  for  green  N2O4  is  -0.061 
-0.063  =  -0.124,  while  for  brown  N2O4  it  is  -0.061  +  0.063  =  0.002.  The 
standard  deviations  of  these  estimates  are  0.04  and  0.016,  respectively. 

Cases  7  and  8  on  the  other  hand  give  estimates  of  -0.046  and  -0.023  for  the 
D  coefficient  for  brown  and  green  N^O^,  respectively.  It  is  possible  that 
the  differences  between  these  cases  could  be  attributed  to  chance.  How¬ 
ever,  it  is  more  probable  that  the  apparent  anomaly  is  due  to  a  modelling 
error.  For  example,  there  may  be  some  cross  effects  due  to  the  initial  tem¬ 
perature;  Cases  7  and  8  were  analyzed  a  chinning  that  there  were  none. 
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This  could  account  for  a  difference  since  T0  =  75  F  for  Case  S  while  Cases  7 
and  8  included  data  with  various  initial  temperatures.  Moreover,  this  argu¬ 
ment  is  not  limited  to  temperature  although  this  is  the  only  parameter  that 
was  held  constant  in  Case  5  and  allowed  to  vary  in  Cases  7  and  8.  Cross  ef¬ 
fects  between  D  and  V  or  D  and  Q/A  could  also  be  present,  causing  biases  to 
be  introduced  since  the  values  for  V  and  Q/A  varied  over  somewhat  different 
ranges  among  the  different  cases. 

When  the  N2O4  is  green,  the  data  indicate  that  the  effect  of  the  filter  his¬ 
tory  parameter,  V,  is  generally  negative  (i.e.,  the  rate  of  flow  decay  de¬ 
creases  as  the  total  amount  of  propellant  which  has  passed  through  the  fil¬ 
ter  is  increased).  In  particular,  Cases  3  and  6  show  this  effect  with  high 
confidence  levels  (0.S9  and  0.84,  respectively).  Case  8  lends  supporting 
evidence  with  a  0.67  level  of  confidence.  Case  S  also  indicates  a  negative 
V  effect  when  no-cross  effects  between  V  and  GB  are  assumed.  However,  it  is 
apparent  from  studying  the  brown  N2O4  runs  that  there  is  a  cross  effect.  In 
fact,  Cases  4  and  7  indicate  that  the  effect  of  parameter  may  be  zero;  both 
compute  V  coefficients  well  within  a/4  of  zero.  Case  1,  on  the  other  hand, 
indicates  a  significant  positive  effect  for  V  when  the  N2O4  is  brown  and 
tested  as  received.  This  case  must  be  considered  cautiously,  since  it  is 
based  on  a  small  number  of  data  points  and  small  errors  in  model  assumptions 
can  have  a  large  effect. 

In  general,  however,  the  investigations  concerning  the  effect  of  V  show  that 
there  is  a  definite  change  in  the  behavior  of  flow  decay  depending  on  whether 
green  or  brown  N^,04  is  being  considered.  Case  5  indicates  that  green  N-,04 
enhances  flow  decay  because  the  GB  coefficient  is  positive.  However,  this 
result  must  be  scrutinized  more  carefully.  It  should  be  recalled  that  all 
the  slopes  were  computed  at  specified  nominal  values  (the  mean  values  for  the 
entire  set  of  test  data)  of  AT,  Q/A,  V,  etc.  This  normalizing  was  done  to 
make  it  possible  to  compare  the  different  cases  directly.  The  problem  is 
that  although  these  are  nominal  values  for  the  entire  data  set,  they  may  be 
atypical  for  a  particular  subset.  The  result  is  that  when  cross  effects 
between  the  variables  are  assumed,  the  main  effects  become  essentially  an 
extrapolation  of  the  data  from  its  own  nominal  value  to  that  of  the  prespeci¬ 
fied  ones.  This  may  not  cause  any  problems  if  the  cross  effects  are  small, 
but  if  they  are  large,  erroneous  results  might  be  computed.  This  may  be  the 
situation  for  Case  5.  A  large  cross  effect  is  computed  between  GB  and  Q/A; 
the  coefficient  of  GB  (Q/A  -  1.5)  is  1.09.  For  this  particular  data  set, 
the  average  value  of  Q/A  is  1.1  and  not  1.5  (the  mean  value  for  the  total 
data  set).  Therefore,  at  the  average  Q/A,  the  effect  of  GB  is  0.41  +  1.09 
(1.1  -  1.5)  =  -0.025.  The  large  coefficient  for  GB  at  Q/A  =  1.5  may  be  the 
result  of  an  extrapolation  error  in  going  from  Q/A  =  1.1  to  Q/A  *  1.5.  The 
general  trend  at  the  nominal  value  of  this  data  set  is  really  negative.  This 
negative  effect  means  that  green  N2O4  experienced  less  flow  decay  than  brown 
N2O4 .  This  conclusion  is  further  substantiated  by  direct  comparisons  between 
sets  of  test  data  for  which  the  only  difference  in  the  nine  basic  independent 
variables  considered  was  between  green  and  red-brown  N2O4. 

In  any  case,  it  is  certainly  apparent  that  the  GB  variable  has  a  large  effect 
on  the  other  variable?.  In  essence,  this  means  that  the  flow  decay  phenom¬ 
enon  behaves  differently  depending  on  whether  green  or  brown  N2O4  is  used. 
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As  another  example,  the  effect  of  the  temperature  parameter  also  appears  to 
be  mixed  depending  upon  whether  or  not  green  or  bi'own  is  used.  For 

green  N2O4,  Case  3  indicates  a  positive  T0  effect  (with  an  estimated  value  of 
0.005  for  the  coefficient  and  0.53  confidence  level  of  a  non- zero  coeffi¬ 
cient).  For  brown  N2O4,  Case  4  predicts  a  negative  effect  with  a  -0.0026 
estimate  and  a  0.47  confidence  level.  Because  of  these  relatively  low  con¬ 
fidence  levels,  these  trends  should  be  considered  with  caution.  Cases  7 
and  8,  on  the  other  hand,  predict  the  same  trends  for  the  TQ  parameter  but 
at  higher  confidence  levels  (-0.019  at  0.999  for  brown  and  +0.013  at 

0.71  for  green  N2O4) . 

The  iron  saturation  parameter  appears  to  have  a  positive  effect  on  flow  de¬ 
cay  for  green  N2O4.  Case  6  substantiates  this  by  estimating  0.40  for  the 
coefficient  of  SAT  and  giving  a  0.87  confidence  level  that  this  coefficient 
is  not  zero.  For  brown  N2O4,  it  is  only  possible  to  make  indirect  compari¬ 
sons  between  as-received  and  doped  propellant.  Regression  Case  1  (brown 
N2O4  tested  as  received)  gives  a  constant  term  of  0.81.  This  term  is  the 
flow  decay  predicted  by  the  regression  curve  at  the  nominal  values  of 
V  =  15,  Q/A  =  1.5,  and  RL  =  0  (i.e.,  representing  an  average  of  the  +1  and 
-1  values  assigned  to  green  and  brown  N2O4,  respectively). 


Case  7  considers  brown,  saturated  N2O4.  The  regression  curve  predicts  a 
value  of  Y  =  1.01,  when  evaluated  at  the  nominal  parameter  values  specified 
by  Case  1  (i.e.,  D  =  5  microns,  AT  =  20  F,  T0  =  75  F,  V  =  15  gal./sq  in., 

Q/A  =1.5  gal./min-sq  in.).  Thus,  there  appears  to  be  a  net  positive  effect 
for  the  saturation  parameter  for  brown  N2O4  as  well  as  green  N2O4 .  Case  4 
might  also  be  used  to  get  a  further  estimate  of  the  saturation  parameter  for 
brown  N2O4 .  The  problem  here  is  that  Case  4  holds  the  diameter  fixed  at  10 
microns,  while  Case  1  holds  it  fixed  at  5  microns,  and  there  is  no  direct  es 
timate  of  the  effect  of  this  diameter  difference  for  these  two  cases. 
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The  regression  curves  indicate  that  the  cooling  rate  parameter  enhances  flow 
decay.  However,  this  effect  is  predicted  only  at  very  low  confidence  levels 
(Case  1  predicts  a  0.003  coefficient  and  a  0.07  confidence  level;  Case  4 
predicts  a  0.12  coefficient  and  a  0.39  confidence  level).  These  low  levels 
are  at  least  partly  due  to  the  fact  that  most  of  the  runs  with  long  cooldowns 
did  not  exhibit  any  significant  flow  decay;  very  few  long  cooldown  runs  with 
significant  rates  of  flow  decay  were  available  for  use  in  the  regression 
analyses. 

The  last  main  independent  variable  that  remains  to  be  considered  is  the  flow 
parameter  Q/A.  The  effect  of  this  parameter  is  generally  positive  (i.e., 
increasing  Q/A  increases  flow  decay).  For  green  N2O4,  Cases  2,  3,  5,  and  6 
predict  this  effect  with  the  very  high  confidence  levels  of  0.93,  0.97, 

0.999,  and  0.99,  respectively.  Case  8  predicts  the  same  effect  at  a  0.60 
level.  For  brown  N2O4,  the  data  are  more  mixed.  Case  1  predicts  a  positive 
effect  for  rapidly  cooled  N2O4 .  This  estimate  is  0.88  with  a  standard  de¬ 
viation  of  0.65.  The  estimate  when  a  long  cooldown  is  used  is  -0.50  with 
a  standard  deviation  of  0.45.  However,  the  results  for  this  case  are  quite 
tenuous  because  of  the  small  number  of  long  cooldown  tests  for  which  flow 
ri n r a v  uac  obc£T"Yscl  ( on ly  Lhyss  points)  ^r.d  Lecnnss  the  1  eT 92  RL  Q/A 
cross  effect. 
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A  deviation  of  just  one  of  the  long  cooldown  points  could  have  a  very  large 
effect  on  the  predicted  equation.  Thus,  modelling  errors,  e.g.,  not  being 
able  to  identify  ail  of  the  independent  variab les,  could  have  a  disastrous 
effect  in  this  particular  case. 

Case  4  predicts  a  negative  Q/A  effect  for  brown -N2O4.  Again,  this  is  only 
at  a  relatively  low  level  of  confidence  (0.55).  The  only  high  confidence 
level  for  the  Q/A  coefficient  for  brown  N2O4  is  given  in  Case  7.  Here  the 
estimate  of  the  coefficient  is  0.62  with  a  level  of  confidence  that  it  is 
not  zero  of  0.99. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

Flow  decay  is  a  very  complex  phenomenon.  There  are  many  independent  variables 
that  affect  the  incidence  and  amount  of  flow  decay.  The  effects  of  these  para¬ 
meters  are  generally  not  simple  or  independent  of  each  other;  they  exhibit  many 
interactions.  In  addition,  there  are  often  threshold  effects  for  flow  decay  (i.e., 
an  identifiable  boundary  between  a  range  of  variables  for  which  no  flow  decay  oc¬ 
curs  and  a  range  for  which  flow  decay  occurs  in  varying  amounts).  These  thres¬ 
holds  are  not  sharp,  and  further  depend  upon  the  interactions  of  the  independent 
variables . 


A  number  of  major  effects  of  parameters  on  flow  decay  are  outlined  in  the  follow¬ 
ing  paragraphs.  The  reader  should  keep  in  mind  the  complex  nature  of  flow  decay, 
and  that  these  are  only  intended  to  be  interim  conclusions.  Extreme  caution  must 
be  taken  in  extrapolating  effects  and  trends  to  different  flow  si  uations  and 
conditions . 

Increasing  the  filter  pore  size  has  a  significant  effect  in  preventing  or  reducing 
flow  decay.  It  appears  that  flow  decay  will  not  occur,  under  the  conditions 
tested  during  this  program,  with  filters  with  nominal  pore  size  at  least  as  large 
as  40  microns.  There  were  a  few  tests  that  might  appear  to  be  exceptions  to  this 
generalization,  however,  they  are  judged  to  be  outliers.  Decreasing  the  filter 
pore  size,  below  its  threshold,  increases  the  rate  of  flow  decay. 

Keeping  the  temperature  drop  (from  the  initial  propellant  temperature)  as  small 
as  possible  is  another  important  factor  in  preventing  or  minimizing  flow  decay. 
However,  some  cases  of  flow  decay  were  observed  with  no  temperature  drop.  There 
is  often  a  non-zero  temperature  drop  threshold  below  which  no  flow  decay  occurs. 

The  amount  or  rate  of  flow  decay  increases  as  the  temperature  drop  is  increased 
above  its  threshold  value.  In  some  ranges  of  conditions,  there  is  approximately 
a  linear  relationship  between  rate  of  flow  decay  and  the  temperature  drop. 

Decreasing  the  local  velocity  through  the  screen  decreases  the  amount  of  flow  de¬ 
cay,  at  least  over  the  range  tested  during  this  program.  The  actual  velocity- 
related  parameter  considered  was  Q/A,  for  which  the  range  of  values  was  about  0.7 
to  4.7  gal/min-sq  in.  There  appear  to  be  strong  interactions  between  Q/A  and 
many  of  the  other  independent  variables.  Use  of  a  filter  Reynolds  Number  was  not 
particularly  advantageous  in  the  data  correlation  efforts. 

Comparisons  between  green  and  red-brown  N?0*  were  not  without  complications. 

There  appeared  to  be  a  number  of  sizable  differences  in  the  effects  of  other  vari¬ 
ables  between  the  two  types  of  propellant,  especially  in  the  effects  of  Q/A,  ini¬ 
tial  temperature,  and  the  history  parameter  V.  However,  direct  comparisons 
between  sets  of  test  data  for  matched  conditions  (i.e.,  sets  for  which  the  only 
difference  in  the  nine  basic  independent  variables  considered  was  between  green 
or  red-brown  N^O^) indicate  that  the  green  N2O,  tested  exhibited  somewhat  lower 
rates  of  flow  decay  than  did  the  red-brown  tested. 

Rapid  temperature  drops  (during  flow)  seem  to  cause  somewhat  more  flow  decay  than 
occurs  with  the  same  temperature  drop  imposed  slowly  over  a  period  of  approximately 
2  days.  However,  this  effect  was  less  important  than  the  effects  of  such 
variables  as  T,  D,  and  Q/A  (all  symbols  are  defined  in  the  Nomenclature  section). 
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The  effect  of  initial  temperature  was  less  consistent  between  differing  types  of 
N->04  when  the  propellants  were  doped  with  iron  pentacarbony  1  to  ensure  saturation. 
With  green  N1O4,  increasing  T0  increased  flow  decay,  while  the  opposite  trend  was 
observed  witfi  red-brown  .\'?04- 

There  was  some  difference  between  the  amount  of  flow  decay  for  the  as-received 
propellants  and  for  the  propellants  doped  with  a  small  amount  (one  saturation 
dose)  of  iron  pentacarbonyl  to  ensure  iron  saturation.  The  differences  seemed  to 
be  regular  and  more  a  matter  of  degree  than  substantive  differences  in  type  of 
flow  decay  behavior.  Therefore,  the  data  for  doped  N204  represent  slightly  con¬ 
servative  results  as  applied  to  the  particular  nitrogen  tctroxide  propellants  that 
were  used  in  the  testing. 

In  addition  to  these  general  effects  of  the  test  conditions,  more  detailed  infor¬ 
mation  is  contained  in  the  body  of  this  report  which  will  be  of  value  to  indivi¬ 
duals  evaluating  the  potential  of  flow  decay  in  operational  nitrogen  tetroxide 
systems.  The  following  portions  of  the  report  should  be  particularly  examined 
for  this  purpose:  Tables  III  through  VII  and  Fig.  4  through  8,  the  text  accom¬ 
panying  these  tables  and  figures,  and  Appendix  A.  For  example,  a  user  could 
select  from  these  test  results  those  which  match  closely  the  conditions  for  a 
system  being  considered;  then,  use  the  results  and  trends  for  those  data  to  help 
in  evaluating  the  flow  decay  potential  for  the  application  of  interest. 


This  program  was  intended  to  be  the  initial  step  in  a  systematic  engineering 
parametric  study  to  establish  the  necessary  engineering  criteria  for  predicting 
flow  decay  in  operational  N204  systems.  Although  many  complexities  were  en¬ 
countered  in  this  initial  step,  it  was  possible  to  develop  experimental  tech¬ 
niques  and  data  analysis  techniques  to  give  statistically  reproducible  flow 

decay  data,  and  to  deduce  from  the  data  the  gross  effects  of  major  parameters. 

It  is  recommended  that  these  efforts  be  continued  to  complete  the  systematic  engi¬ 
neering  parametric  study.  Subsequent  efforts  would  build  upon  the  results  of 
this  program  and  culminate  in  a  design  guide  for  predicting  the  occurrence  or 
absence  of  flow  decay  and  establishing  design  and  system  management  criteria  to 
avoid  this  problem. 

The  tests  during  this  program  were  intended  to  investigate  gross  effects  of  those 
parameters  that  were  expected  to  have  the  greatest  influence  on  ^^4  flow  decay. 
The  recommended  future  tests  woulu  extend  the  work  to  delimit  more  precisely  the 
ranges  of  these  operating  and  design  parameters  under  which  flow  decay  will  and 
will  not  occur,  determine  the  effects  of  additional  parameters  (e.g.,  concentra¬ 
tions  of  NO,  H20,  and  other  impurities  in  the  N2O4,  plus  other  types  of  flow 
constrictions),  and  investigate  more  thoroughly  the  interactions  between  various 
parameters.  Completion  of  this  systematic  engineering  study  would  form  the  basis 
for  engineering  control  of  this  potential  system  failure  mode  by  defining  required 
system  design  criteria  and/or  system  management  concepts  to  avoid  flow  decay. 
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APPENDIX;  COMPUTER  PRINTOUTS  OF  DERIVED  RESULTS 


The  computer  output  given  in  this  appendix  lists  for  each  test  run  values  of 
several  derived  parameters  for  each  of  the  possible  pairings  of  the  data 
points  taken  from  the  recorder  traces.  The  particular  data  points  are  in¬ 
dicated  by  the  two  time  values  shown  in  the  second  column.  The  third  and 
fourth  columns  list  values  of  X  and  Y  (defined,  with  other  symbols,  in  the 
text  and  in  the  Nomenclature  section).  The  fifth  column  lists  values  of  W 
for  each  of  the  time  values  given  in  the  second  column.  The  abbreviated 
labels  near  the  right  side  of  the  table  give,  in  this  order:  (1)  type  of 
N2O4  ("B"  for  red-brown  and  "G"  for  green),  (2)  saturated  ("S")  with  iron  or 
tested  as  received  (left  blank) ,  (3)  cooled  rapidly  ("R")  or  over  a  long 
time  ("L"),  (4)  the  filter  size  (normal  pore  size  in  microns,  followed  by 
"F"  or  "H"  to  distinguish  different  filter  areas;  generally,  the  filters 
marked  MH"  had  roughly  one-half  the  area  of  filters  marked  "F",  although 
actual  areas  were  used  in  all  calculations),  (5)  the  initial  temperature, 
and  (6)  the  temperature  drop.  For  Runs  1-38  and  160-165  calibration  tests 
were  not  made  for  the  clean  filters.  Consequently,  wl ,  w2 ,  2,  and  1/L  are 
not  defined.  This  fact  is  indicated  by  the  abbreviation  "N/A"  in  the  appro¬ 
priate  locations  in  the  tables. 
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30  (J«6/  3*4  -0.14B  -0.053  -0.40/ -0-41  -0*002  3  909900  2.Q  BS  R  5F  75  0 

30  0*6/ 6*2  0-400  0.073  -0.40  /  0*00  0.072  6  1513  1.9 

30  3*4/  6«2  Q.556  0.19B  -0-41/  04)0  0*147  3  552  1*9 

40  043/3.4  14533  0*513  -0450/  0*79  0-455  3  224  4.1  BS  R  10F  75  0 

40  043/042  1.075  0-190  -0459/  0-10  0.140  6  571  44) 

40  3*4/  6«2  -0.261  -04)93  0.79/  0-40  -0.130  3  999993  443 

41  0.67  2*6  0.300  0*150  0-37/  0.74  0.195  2  320  2.2  BS  R  5F  75  0 

41  0.6/  6.6  0*636  J.106  0.37/  1-12  0.126  7  1151  2.1 

41  2*6/  6*6  0-337  0-034  0.74/  1.12  0.096  5  1443  2.1 

42  1*4/  4*4  0*522  0*174  0*00/  0*36  0*120  4  726  4*4  BS  A  10F  7b  O 

4<i  1*4/  7*0  0*013  0*163  0-00/  1*00  0*195  7  770  4*4 

4i>  4.4/  7.0  0*392  0*151  0*36/  1*09  0.231  3  325  4.4 

43  2*4/  7.2  1.301  0.200  2.50/  4.06  0>306  6  422  1.0  BS  K  6F  7b  21 

43  2-4/13.6  3-370  0.316  2»50/  6.©?  0*365  15  340  1*0 

43  7.2/13.6  2.523  0-394  4.0iV  6.67  0.407  7  303  1-3 

44  0.2/  6*0  1*736  0.290  1.11/  2.03  0.310  7  406  3.7  j>S  it  10F  75  21 

44  0-2/14.0  3-510  0-254  1.11/  4*51  0*246  16  471  3.7 

44  3-0/14-0  1*005  0*22o  2*9o/  4*ol  0*104  0  525  3*  / 

45  1.0/  7.2  4.739  0.772  9.02/  14.00  O-oi?  7  159  1-6  BS  K  5F  74  35 

45  1.0/14.0  13-310  0.964  9.02/  21.55  0.350  15  123  1*5 

45  7.2/14.0  0-050  1.170  14-00/  21-55  0.004  0  00  1*5 

46  1.0/  <3.0  2.590  0-470  4.02/  7.20  0-455  6  240  3-3  33  R  10F  74  55 

46  1.0/12.2  6-666  0.594  4-02/  11.03  0.5S7  1.2  196  3.2 

46  1.0/10.0  11.031  0.702  4.02/  16.30  0*660  10  162  3.1 

46  6-0/12.2  4.369  0.705  7.20/  11-23  Q.6E0  7  133  5.2 

46  0.0/10*0  0.774  0*015  7-20/  16*30  0*70  13  130  3*1 

4o  12.2/14-0  5-652  0*974  11*29/  15-30  0-665  6  113  3.0 
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